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We argue that the typical energy density of a light scalar field should not be less than H 4 in the 
inflationary Universe. This requirement implies that the non-Gaussianity parameter /jvl is typically 
bounded by the tensor-scalar ratio r from above, namely /atl i; 5f8 • r*. If jnl ~ 10 , inflation 
occurred around the GUT scale. 
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I. INTRODUCTION 

Inflation model 

0,0,1 provides an elegant mechanism 
to solve the horizon, flatness and primordial monopole 
problem due to a quasi-exponential expansion of the 
universe before the hot big bang. The temperature 
anisotropics in cosmic microwave background radiation 
(CMBR) and the large-scale structure of the Universe 
are seeded by theprimordial quantum fluctuations dur- 
ing inflation 0, (5[. Since the density perturbation is 
roughly 1CU 5 , it is good enough to apply the linear 
perturbation theory to calculate the quantum fluctua- 
tions during inflationary epoch. Within this approach, 
the Fourier components of fluctuations are uncorrelated 
and their distribution is Gaussian. That is why the 
non-Gaussianity from the simplest inflation models is 
very small {\/nl\< IV For useful discussions on non- 
Gaussianity see [& 0, B i, [M G3 , 

and for a nice review 

see Q3. 

The non-Gaussian perturbation is governed by the n- 
point correlation function of the curvature perturbation 

n 

($(k 1 )$(k 2 )-$(k n )) = (2 7 r) 3 5 3 (^k i ) J F n (fc 1 ,A: 2 ,...,fc„), 

i=l 

(1) 

where $(k) is the Fourier mode of Bardeen's curva- 
ture perturbation. The leading non-Gaussian features 
are known as the bispectrum (three-point function) and 
trispectrum (four-point function), with their sizes con- 
ventionally denoted as $nl and t^l respectively. The 
non-linearity parameter /jvl defined in [l3T ] takes the 
form 



*(x) = * L (x) + /jvl[*£(x) - ($ 2 L (x))], 



(2) 



where denotes the linear Gaussian part of the per- 
turbation in real space. The non-Gaussianity parameter 
/jvl characterizes the amplitude of the primordial non- 
Gaussian perturbations. The value of <&z, is roughly 10~ 5 . 
If /jvl = 10 2 , the distribution of <E> is still consistent with 
a Gaussian distribution to 0.1%. It is not easy to detect a 



small non-Gaussianity in the experiments. The shape of 
non-Gaussianity is also very important. If F(ki, hz, ks) 
is large for the configurations in which k% <C fc 2 , k% , it 
is called local, "squeezed" type; if F(ki, k2, k^) is large 
for the configurations in which k\ ~ fc 2 ~ &3, it is called 
non-local, "equilateral" type. 

Recently a large primordial local-type non-Gaussianity 
was reported to be marginally detected from the data of 
WMAP3 in [3: 



27 < f l N C L al < 147 



(3) 



at 95% C.L., with a central value of Jml = 87. See also 
(l5| . In [l6[ WMAP group recently reported their new 
result as 



-9</#£ a/ <m- 



(4) 



If a large local-type non-Gaussianity is confirmed by the 
future cosmological observations and data analysis at 
high confidence level, it strongly implies that the physics 
of the early Universe is more complicated than the simple 
single-field slow-roll inflation. 

There are several mechanisms to generate a large local- 
type non-Gaussianity: 

(1) a non- linear relation between inflaton and curvature 
perturbations 0; 
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(2) curvaton scenario 

(3) Ekpyrotic model 
In this paper we focus on curvaton scenario. As we know, 
there might be many light scalar fields in supersymmetric 
theory or string theory. These light scalar fields which are 
subdominant during inflation are called curvaton. The 
energy density during inflation is still dominated by the 
the potential of inflaton. In the usual inflation model 
the fluctuations of inflaton dominate the curvature per- 
turbations. The energy density and the perturbations 
caused by these light scalar fields can be ignored during 
inflation. However it is possible that the fluctuation of 
curvaton becomes relevant and causes a large local-type 
non-Gaussianity when its energy density is a significant 
fraction of the total energy after the end of inflation. 
In curvaton scenario the perturbations from the infla- 
ton field are considered to be negligible. A large local- 
type non-Gaussianity may shed light on these light scalar 
fields. 
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The non-Gaussianity produced by curvaton is in in- 
verse proportion to the fraction of curvaton energy den- 
sity in the energy budget at the epoch of curvaton de- 
cay. The non-Gaussianity increases as the curvaton en- 
ergy density decreases. However the energy density for a 
homogeneous scalar field in the quasi-de Sitter space is 
typically is given by Hf. Therefore the non-Gaussianity 
has an upper bound 



Inl < 518 • 



(5) 



where r is the tensor-scalar ratio. Present bound on the 
tensor-scalar ration |16] is r < 0.20 and then f^L ^ 346. 

Our paper is organized as follows. We review curva- 
ton scenario and the curvature perturbation generated 
by curvaton in Sec. 2. In Sec. 3, we argue that the 
energy density of curvaton is typically estimated as Hf 
during inflation and then the non-Gaussianity parameter 
/jvl is bounded by the tensor-scalar ratio from above. In 
Sec. 4, we point out the challenge for getting a red-tilted 
primordial power spectrum in curvaton scenario. Section 
5 contains some discussions including how to distinguish 
the curvaton scenario from Ekpyrotic scenario. 



II. NON-GAUSSIANITY IN CURVATON 
SCENARIO 

In curvaton scenario the dynamics of inflation in the 
early Universe is dominated by inflaton and the en- 
ergy density of curvaton is negligible during inflation. 
The lagrangian of a curvaton field a in an unperturbed 
Friedmann-Robertson- Walker (FRW) spacetime is given 
by 



1 



1 



-(V p a) 2 -V(a), 



(6) 



where the subscript "p" denotes the derivative with re- 
spect to the physical coordinate. For simplicity, we con- 
sider a potential of curvaton as follows 



V(a) = \m 2 a 2 . 



(7) 



The curvaton field is supposed to be an almost free field 
with a small effective mass compared to the Hubble scale 
H during inflation. Effectively we define a parameter 



1 d 2 V(a) 
3H 2 da 2 



(8) 



which is much smaller than 1. The equation of motion 
for a homogeneous curvaton field takes the form 



a + 3Ha 



-m 2 a. 



(9) 



Since m <C H, the friction term 3Ha is dominant and the 
value of curvaton is roughly a constant at the inflationary 
epoch. It is denotes as a = er*. 



The amplitude of the quantum fluctuation of curvaton 
field in a quasi-de Sitter space is given by 



da — — . 

2tt 



(10) 



The spectrum of the fractional perturbations caused by 
curvaton fijl is 



p-2 

8a I 'a 



6a 
a 



2tt a* 



(11) 



where * denotes the epoch of horizon exits fc* = a*H* 
during inflation. The quantum fluctuation of curvaton 
field is frozen at horizon exit to a classical perturbation 
with a flat spectrum. Since the curvaton energy density 
is subdominant in this epoch, its fluctuations are initially 
taken as isocurvature/entropy fluctuations. But the cur- 
vature perturbation is sourced on large scales [12J. 

After the end of inflation the inflaton energy density is 
converted into radiations and then the Hubble parameter 



decreases as 



The curvaton field will remain approx- 



imately constant <r* until H ~ m. At this epoch the cur- 
vaton starts to oscillate harmonically about a = 0. Dur- 
ing the stage of oscillating the curvaton energy density 



goes like p a 



which increases with respect to the 



energy density of radiation pn tx a" 4 . When the Hubble 
parameter goes to the same order of the curvaton decay 
rate T, the curvaton energy is converted into radiations. 
Finally, before primordial nucleosynthesis, the curvaton 
field is supposed to completely decay into thermalized ra- 
diation, thus the perturbations in the curvaton field are 
converted into curvature perturbations and become the 
final adiabatic perturbations which seed the matter and 
radiation density fluctuations observed in the Universe. 
We illustrate the evolution of one curvaton perturbation 
mode and Hubble radius in Fig. 1. 




oscillates 



Curvaton 
decays 



FIG. 1: Evolution of the perturbation mode and the Hubble 
radius. 

After the curvaton decay the total curvature pertur- 
bation will remain constant on superhorizon scales at 
a value which is fixed at the epoch of curvaton decay. 
The sudden-decay approximation is exact if the curva- 
ton completely dominates the energy density before it 
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decays. Going beyond the sudden-decay approximation, 
a numerical calculation in [271 ] implies that the spectrum 
of the curvature perturbations caused by curvaton is pro- 
portional to the fraction of curvaton energy density in the 
energy budget at the epoch of curvaton decay 



(-) 

\PtotJ , 



(12) 



A precise form of the spectrum in curvaton sceanrio is 
given in [2l| as follows 

p? =~n aD p? , =—n aD —. (13) 

In curvaton scenario the curvature perturbation is dom- 
inated by curvaton perturbation and then the spectral 
index in the curvaton model takes the form 

n s = 1 + 2r/ CTCT - 2e#, (14) 

where ch is the slow-roll parameter which is defined as 

/ / 



(15) 



In [2l|,[28[ the non-Gaussianity parameters corresponding 
to bispectrum 1 and trispectrum are respectively given by 



and 



INL 



TNL 



4n„ 



36 
25 



NL- 



(16) 



(17) 



In sudden-decay approximation, Eq. (|16p should be re- 
placed by 5/(30 o - i _d). A large non-Gaussianity is ob- 
tained if fl a .D <C 1. In this case the Hubble parameter is 
always dominated by the radiation energy density before 
curvaton decays. Assume the scale factor is a{t ) = 1 at 
the moment when curvaton start to oscillate. The energy 
density of curvaton and radiation are p a {t ) = ^m 2 a 2 
and pn{t ) — 3M 2 m 2 respectively. Since the radia- 
tion energy density goes like a -4 , pR(tn) — 3M 2 r 2 ~ 
P_R(*o)a 4 (^r>) an d thus a(to) ~ (m/T)i. Therefore 



.D 



D 



err [m\2 



p R {t D ) %M1 



(18) 



This result is also given in [l9|. Combing ([16]) with (fjg]) , 
we get 



15M 2 

JNL = -7. T — 



(19) 



1 A precise form in [T3l is /j\ 



Since T < m, a large non-Gaussianity is obtained only 
when er» -C M p . It is reasonable to require that the VEV 
of curvaton is less than Planck scale and a large Jnl is 
expected. 

The WMAP normalization is 



P c = 2.457 x 10" 



(20) 

Using eq. (TIB"]) and (flB)) . the value of curvaton during 
inflation is related to the Hubble parameter by 



cr* = 2.68 x 10-* 



NL 



(21) 



The non-Gaussianity parameter /jvl should be smaller 
than 2.68 x 10 3 ; otherwise, the quantum fluctuation of 
curvaton is greater than its VEV and the previous semi- 
classical description breaks down. Substituting (|21jl into 
(Pj). we find 



NL 



10 ! 



H 2 /my 



(22) 



The limit of the primordial gravitational wave pertur- 
bation implies that < 10~ 4 M p , and then the mass 
of curvaton must be much larger than its decay rate for 
getting a large non-Gaussianity. 

Tensor/gravitational wave perturbation encodes a very 
important information about inflation. Its amplitude 
only depends on inflation scale 



Pt 



Hi/Ml 

7T 2 /2 



(23) 



For convenience, we define a new parameter named the 
tensor-scalar ratio r as 



Pi 
P t 



(24) 



If the curvature perturbation is dominated by inflaton 
field we have r = 16e# for slow-roll inflation. In cur- 
vaton scenario the density perturbation is dominated by 
curvaton and thus 



r < 16en- 



(25) 



According to (|23|) , the Hubble parameter during inflation 
is related to the tensor-scalar ratio by 



/ / 



Now eq. 



= 10- 4 
becomes 



r^Mp = 2.44 x 10 1 



■^GeV. 



Jnl = 10" 



r 



(26) 



(27) 



WMAP 5 gives an upper bound on the tensor-scalar ratio 
r < 0.20 PJ] and we need m/T > 2.5 x 10 9 for f NL = 10 2 . 
The bound on the tensor perturbation also provides an 
upper bound on the inflation scale 

H* < 1.09 x 10 14 GeV. (28) 

We will see that the non-Gaussianity in curvaton scenario 
will give a lower bound on H* in the next section. 
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III. UPPER BOUND ON NON-GAUSSIANITY 
f NL IN CURVATON SCENARIO 

The curvaton mass and its decay rate are fixed if the 
field theory for the whole system is given. Since the 
curvaton does not move during inflation, we cannot use 
its dynamics to determine the value of a*. According 
to Eq. (flH|) , /jvl is large if cr* is much smaller than the 
Planck scale. In the literatures, tr» is not fixed by the the- 
ory, but rather represents an additional free parameter, 
and then the non-Gaussianity parameter /jvl can be ar- 
bitrarily large. This treatment is reliable at the classical 
level. However we will see that the quantum fluctuations 
of curvaton will significantly affect the value of curvaton 
during inflation and an upper bound on the typical value 
of fNL is obtained for the curvaton model. 

In the inflationary Universe the quantum fluctuations 
of the curvaton field freeze out with amplitude H/ (2?r) 
at the wavelength H . In [29j a gradient energy density 
due to these fluctuations is estimated 



(V p a) 2 



(29) 



It is also obtained through calculating the contribution 
of these quantum fluctuations to the average value of 
the energy momentum tensor. The gradient energy den- 
sity characterizes the inhomogeneity of the curvaton and 
should be smaller than its homogeneous energy density 
^m 2 cr 2 . Requiring m 2 a 2 > H% yields a lower bound on 
the VEV of curvaton field, namely 



> 



(30) 



Since m < H„, er» > H* which is the condition for the 
validity of the semi-classical description. 

Secondly lets take into account the horizon tempera- 
ture £/* of a quasi-de Sitter space. Since the curvaton 
mass is much smaller than iJ*, the curvaton "particles" 
can be taken as the relativistic particles. If the initial 
temperature of curvaton particles is lower than if*, the 
quantum fluctuations of the background will automati- 
cally heat their temperature up to H*. So the energy 
density of curvaton field should not be smaller than H%. 

Here we also give the third argument. In curvaton 
model the curvaton mass is assumed to be much smaller 
than the Hubble parameter during inflation, which means 
the Compton wavelength is large compared to the cur- 
vature radius of the de Sitter space -ff" 1 . So the grav- 
itational effects may play a crucial role on the behavior 
of the curvaton field. The well-known Bunch-Davies ex- 
pression [10, HH, H3 for the average square value of a light 
scalar field (m *C i/*) in a quasi-de Sitter space takes the 
form 



3Ht 
87r 2 m 2 



(31) 



An intuitive understanding was given in [33j |. According 
to the long-wave quantum fluctuation of a light scalar 



field (m <C H) in inflationary universe, the behavior of 
such a light scalar field can be taken as a random walk 



^ = 4^ 



(32) 



On the other hand, a massive scalar field cannot grow 
up to arbitrary large vacuum expectation value because 
it has a potential. The long wavelength modes of the 
light scalar field are in the slow-roll regime and obey the 
slow-roll equation of motion, i.e. 



dt 



dV(a) 
da 



-m 2 a. 



(33) 



Combining these two considerations, in [33| Linde and 
Mukhanov proposed 



d(a 2 
dt 



4tt 2 



2m 2 

w 



(34) 



We see that (a 2 ) stabilizes at the point of (a 2 ) 
which is the same as Eg. pip . This vacuum expectation 
value of curvaton mainly comes from the perturbation 

H 2 _i 

mode with wavelength — fif, which stretched outside 
the horizon at the number of e-folds H 2 / m 2 . The infla- 
tion models with small total number of e-folds are arti- 
ficial and the long stage of inflation is expected generi- 
cally [iH, [36[ . Since the wavelength is much larger than 
the Hubble horizon, this fluctuation mode is frozen to be 
a classical one and provides a non-zero classical configu- 
ration for curvaton field. The typical value of curvaton 
field in such a background is estimated as 



Hi 

m 



(35) 



Similar estimation is also adopted in the literatures, for 
example (37| . 

The value of er» is related to the Hubble parameter 
and /jvl through the WMAP normalization by eq. (|2ip . 
Eq.([3SD reads 



m 



NL 



2.68 x 10 3 ' 



(36) 



the 



Because the mass of curvaton is smaller than H r 
non-Gaussianity parameter is less than 2.68 x 10 3 . 

The curvaton should decay before neutrino decoupling; 
otherwise, the curvature perturbations may be accompa- 
nied by a significant isocurvature neutrino perturbation. 
The temperature of the universe at the moment of neu- 
trino decoupling is roughly T n d — 1 MeV. So the curvaton 
decay rate is bounded by the Hubble parameter at the 
moment of neutrino decoupling, namely 



r>r 



T 2 



4.1 x 10 _25 GeV 



(37) 
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Considering |22J and m < H*, we find 



Inl < 10 6 



Hi 



Similarly combining 



H* 



f. 



NL 



< 



2.71 x 10 7 GeV 
and (j3"6")l , we have 

5 

m \ 6 



(38) 



4.9 x 10 4 GeV 



(39) 



For a large non-Gaussianity Jml — 10 2 , the bounds on 
the Hubble parameter and the mass of curvaton become 

> 1.71 x 10 8 GeV, and m > 1.23 x 10 7 GcV. (40) 

The above constraints on the Hubble parameter and cur- 
vaton mass are not restricted. 

In [H, [39[ the authors also suggested that the curvaton 
decay rate is at least of order Gm 3 ~ m 3 /M 2 correspond- 
ing to gravitational-strength interactions. So we have 



r>r ~ — 

~ 9 ~ Mr 



Now eq. (|22j) reads 



toM b 



Including (|36p . we find 



fNL < 



3.4 x 10 5 GeV 



(41) 



(42) 



(43) 



If /jvx = 10 2 the curvaton mass should be larger than 
3.4 x 10 11 GeV. Combing (SB), gU) and flSBJ), we find 
the typical values of /jvl and tjvl are bounded by the 
tensor-scalar ratio r from above: 



Or equivalently, 



f NL < 518 Ti, 
tjvl < 3.86-10 5 -ri 

ff» > 9 • 10 8 • f^r L GeV. 



(44) 
(45) 

(46) 



Present bound on the tensor-scalar ratio is r < 0.20 and 
then f NL < 346. If f NL = 10 2 , r > 1.4 X 10~ 3 , ff* > 
9 x 10 12 GeV and inflation scale V" 1 / 4 > 6.2 x 10 15 GeV 
which implies that inflation happened around GUT scale. 
The inequality (|44|) is saturated only when Eq.([4"Tj) are 
saturated. So the decay rate of curvaton cannot be much 
larger than gravitational-strength decay; otherwise the 
non-Gaussianity will be small. 

The present and next generation of experiments, such 
as WMAP and PLANCK, will increase the accuracy to 
about ~ 6 and \tnl\ ~ 560. If there is really a 

large local-type non-Gaussianity, it will be detected at 
high confidence level in the near future, and eq. (|4~4"|) 
provides a restricted constraint on the curvaton scenario. 



IV. SPECTRAL INDEX OF THE PRIMORDIAL 
POWER SPECTRUM 

Spectral index of the primordial power spectrum is an- 
other important quantity for characterizing the physics 
in the early Universe. WMAP 3 [l6[ prefers a red tilted 
primordial power spectrum with n s — 0.9601q'oi3 for 
ACDM model, and n s > 1 is disfavored even when grav- 
itational waves (r < 0.2 at 95% CL) are included. 

The spectral index in the curvaton scenario is given by 

2 m 2 

n s -l = 2i laa -2t H = -—-2e H . (47) 

According to eq. ([56]) . we see that r] aa is bounded by the 
non-Gaussianity parameter /jvl from below, namely 



Vo 



4.64 • 10" 8 • / 



2 

NL- 



(48) 



In curvaton scenario the curvature perturbation is dom- 
inated by curvaton perturbation and then the slow-roll 
parameter e# has a lower bound (|25[) which says 



e H > — > 8.68- 10 
16 



-13 



Inl- 



(49) 



If Inl = 10 2 , 2r] aa - 9.28 • 10~ 4 and 2e H > 1.74 ■ 10~ 4 . 
We see that the non-Gaussianity cannot give us a useful 
constraint on the spectral index. 

Usually we assume m <§; H * and then r] acr ~ 0. To get 
the observed spectral index we need ch — 0.02 in cur- 
vaton scenario. Among so many inflation models, large- 
field models, such as chaotic inflation, give roughly the 
suitable value of en- Large field means the VEV of infla- 
ton is larger than the Planck scale. However in El| 
we argued that the VEV of a scalar field should be less 
than the Planck scale in a consistent low-energy effec- 
tive field theory coupled to gravity. On the other hand, 
the Lyth bound [42j for single-field slow-roll inflation is 
roughly given by 



(50) 



where N e is the number of e-fold before the end of infla- 
tion. Requiring \A<j>\/M p < 1 yields 



e H < 



2(AN e ) 



(51) 



For AiV e = 50, en < 2 x 10~ 4 . So it is reasonable to 
expect that a spectral index indistinguishable from 1 is 
obtained in curvaton scenario. This is also pointed out 
in [43| . In inflation model (n s — 1) is mainly due to a 

negative r)^ = M^^M^. For the detail see 0. 

A possible way to avoid the above bound is to consider 
multi-field inflation model. The simplest one is assisted 
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inflation. Since there is a unique attractor behavio 
— — the generalized Lyth bound reads 



Z>2 = ■■■ = Pn, 



A G 



= V2^AN e , i = l,2,...,N, 



ir 0i = 



(52) 



where Aq = M p /y N is the gravity scale for N species 
0, [45[. In (44| we gave several examples to support 
that the variation of each inflaton should be less than 
Aq in the assisted inflation model. If so, the multi-field 
inflation cannot help us to release the constraint on the 
slow- roll parameter en- 

Another possibility is chain inflation in string land- 
scape [!(| In this scenario, the universe tunneled 
rapidly through a series of metastable vacua with dif- 
ferent vacuum energies. Since the total energy density 
is dominated by the vacuum energies, we get an infla- 
tionary universe. Chain inflation is not really a slow-roll 
inflation model and it does not suffer from the above 
constraints. In [47| we have e# = j^-- For N e = 50, 
(H = 0.0067, n s = 0.987. Many light scalar fields are ex- 
pected to emerge in four-dimensional effective field the- 
ory from string theory. Some of them can be taken as 
curvatons. The bound on the non-Gaussianity parame- 
ter is Jnl ^ 296. Chain inflation might be generic in 
string landscape. So the experiment data can be nicely 
explained in string landscape. 

Because the curvature perturbation is dominated by 
curvaton perturbation in curvaton scenario, we have r < 
1Q(h and 



Inl < 871 



\AN e M p J 



(53) 



here we take (gU) and ([50]) into account. If \A(j)\/M p < 1 
for AN e — 50, Jnl ^ 123 which is roughly saturated the 
experiment bound. 



V. DISCUSSIONS 

An unambiguous detection of Jnl > 10 will rule out 
most of the existing inflation models. The uncertainty of 
fx c L al from WMAP will shrink to be 42 for 8 years data, 
and 38 for 12 years data The accuracy of PLANCK 
will be roughly 6. If the local- type non-Gaussianity is of 
order 10 2 , it will be detected at high confidence level in 
the near future. 



In this paper we investigate the curvaton scenario in 
detail and find that the inflation scale is roughly at GUT 
scale in order to get a large local-type non-Gaussianity. 
Ekpyrotic model can also provide a large local-type non- 
Gaussianity. However unlike the slightly red-tilted grav- 
itational wave spectrum in the inflation/curvaton model, 
the gravitational wave spectrum in Ekpyrotic model is 
strongly blue and then the amplitude is exponentially 
suppressed on all observable scales 48]. Gravitational 
wave perturbation can be used to distinguish curvaton 
scenario from Ekpyrotic model. 

Here we also want to clarify two points in this paper. 
One is that maybe we miss some order-one coefficients 
in ([35]) and (j4TJ), and then the bound in ([44]) should be 
modified to be a little looser or more stringent. But the 
order of magnitude of our result can be trusted. The 
other is that the bound in (|44]) depends on the WMAP 
normalization P^.wmap = 2.457 x 10~ 9 . If we let the 
normalization of the density perturbation free, we find 



f N L < I' 1 



518-r* 



(54) 



and eq. ([5"3"]) becomes 



fNL 



<7 _1 -871- 



AJV e M p J 



(55) 



where 7 = (P^ /P$,wmap)*- Larger the amplitude of 
the density perturbation, smaller the non-Gaussianity. If 
we also consider \A<f>\/M p < 1 for AN e = 50, /jvx i5 
123 • 7 _1 . The density perturbation cannot be larger 
than 10 times of its observed value in our Universe if the 
non-Gaussianity parameter is not smaller than 10. 

The equilateral-type non-Gaussianity has not been de- 
tected. Many models [5(1 [EH, [H, [53[ were sug- 
gested to generate large equilateral-type non-Gaussianity 
as well. Other mechanisms concer ning the large non- 
Gaussianity are discussed in [H, [HI, [Ha, H3, [H, H^, [6(| • 
Anyway, non-Gaussianity will be an important issue for 
cosmology and string theory. 
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